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Abstract 

We use the exact scattering description of the scahng Ashkin- Teller model in two 
dimensions to compute the two-particle form factors of the relevant operators. These 
provide an approximation for the correlation functions whose accuracy is tested 
against exact sum rules. 



1. The computation of correlation functions is a central and notoriously difficult prob- 
lem of quantum field theory and statistical mechanics. The integrable models of two- 
dimensional quantum field theory provide an unique framework where such a problem 
can be studied in a non-perturbative way. Although even in this context the determi- 
nation of exact analytic expressions for the correlations functions remains a challenging 
task, it is remarkable that an extremely accurate evaluation of the correlators is nowdays 
possible and involves a relatively small amount of technical work. The method is quite 
general and relies on the solution of the integrable model in terms of an exact scattering 
theory [Q. The natural way to express correlators is then as spectral series over a com- 
plete set of intermediate asymptotic states. The computation of the matrix elements of 
the operators (form factors) entering the spectral series is the crucial intermediate step 
needed to make contact with the space of local operators starting from the purely on-shell 
solution of the model [|], |^. From the physical point of view, a particularly interesting 
point in the problem of the computation of form factors lies in the identification of the 
solution of the form factor equations corresponding to a specific operator. In Refs. |^, 
a set of operator-dependent constraints for the form factors was identified which allows 
the solution of this identification problem. The need for approximation in this program 
for the computation of correlators arises only from our present inability to resum exactly 
the spectral series, what forces us to rely on partial sums. 

In this note we briefiy illustrate how this set of ideas applies to the scaling limit of 
the Ashkin- Teller (AT) model. The richness of the spectrum of relevant operators, the 
presence of different locality sectors and the relation with the Sine-Gordon model are 
among the features which make this example interesting from the theoretical point of 
view. 

The two-dimensional AT model can be formulated in terms of two planar Ising 
models coupled by a local four-spin interaction. The lattice Hamiltonian reads 

Hat = E + ^>2) + i^«ar^2] , (1) 

(m,Ti) 

where cr™2 = =tl are the two Ising spins at the site m, the sum is over nearest neighbours 
and the same coupling J has been chosen for the Ising models (isotropic case). The 
model is known to exhibit a line of second order phase transition in the space of the 
lattice couplings J and K along which the critical exponents vary continously |^. In 
this note we will consider the scaling region around this critical line. In such a region 
the correlation length is much larger than the lattice spacing and the model admits a 
continous description in terms of the Hamiltonian 

Hat = n'i + nl + T j (fx{ei{x)+e2{x))+ p j (fxei{x)e2{x) , (2) 
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where and denote the fixed point Hamiltonian and the energy density operator 
of the i-th Ising model, respectively. The meaning of the last expression is that the 
scaling limit of the AT model can be regarded as a conformal field theory with twice the 
central charge of the Ising model (i.e. c = 1/2 + 1/2 = 1) perturbed by the two operators 
S = El + 62 and 6162. 

Consider at first the two critical, non-interacting Ising models, namely the Hamiltonian 
with T = p = 0. The operator 6162 has twice the scaling dimension of the Ising energy 
density, namely x^^ej = 2. Then it is marginal and taking p 7^ in (H) does not spoil 
criticality. This means that the coupling p parameterises the critical line of the AT model. 
A critical line with central charge c = 1 can be described in terms of a free massless boson 
ip (Gaussian model). The correspondence between the critical AT model and the Gaussian 
model was established in Ref. (see also |T^) where, in particular, a series of operator 
identifications was obtained. These identifications are summarised in the first two columns 
of the Table. 

We denote by cri{x) and Pi{x) {i = 1,2) the order (spin) and disorder operators of 
the two Ising models. It is a characteristic feature of the AT model that these operators 
preserve along the whole critical line the basic properties they have in the pure Ising 
model 0]. In particular, they retain the same scaling dimension Xo- = = 1/8, and 
produce the neutral fermions ipi under operator product expansion: at x pi ipi. 

The parameter /3, equivalent to p, identifies the different points along the critical line 
in the Gaussian model. The field (p, dual to 99, is defined by idaP = Eapdisip, with a = 1, 2 
labelling the two directions on the plane. Within the normalisation of the bosonic field ip 
we adopt in this paper, the vertex operators e*""^ and e^°"^ have scaling dimension a'^/An, 
what completes the third column of the Table. The point p = at which the two Ising 
models are decoupled corresponds to P = -y/ivr. It was argued in |^ that the AT critical 
line corresponds to the range \/2ti < (3 < -\/87r. 

When r 7^ in the Hamiltonian (0) the system is moved away from criticality and 
develops a finite correlation length. In view of the identification £ ~ cos j3p, we see that 
the bosonic equivalent of the resulting massive theory is the Sine-Gordon model defined 
by the action 

AsG = J d^x i^df^pd^p + pcos(3p) , (3) 

where p is a dimensional parameter proportional to r. The equivalence of the models 
(^ and (^ (for the range of P specified above) immediately enables us to extend to the 
scaling AT model the two basic results known for Sine-Gordon, i.e. its integrability and 
its exact scattering description |]T|. 

The exact particle spectrum of the Sine-Gordon model contains two elementary excita- 
tions of mass m, the soliton and the antisoliton A_. If from the bosonic point of view 
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they are topological excitations, they can also be interpreted as the particles associated 
to the charged elementary fermion of the equivalent massive Thirring model |Tl|. If we 
write A± = Ai ± iA2, then the neutral excitations Ai and A2 are naturally associated to 
the neutral fermions ipi and ijj2 of the AT model. 

Due to integrability, the scattering of the particles Ai and A2 is completely elastic and 
factorised. It is entirely specified by the Faddev-Zamolodchikov algebraQ 

AMM02)= E a'[j{e,-92)Ai92)Aki9,) . (4) 

k,l=l,2 

The non-zero scattering amplitudes are given by |I|| 

2all{e) = 2all{e) = s{e) + s4e) , 

2alji9) = ^^11(0) =S+{e)- 3(9) , 

2aim = 2allie) = sie) + s4e) , 

2ajl{e) = 2all{e) = s{e)-s4e) , (5) 

where 

sinh ^(6 + m) cosh ^A6 + m) 




<x> rfa; sinh f (1 - f ) Ox 

—I I T sm — 

X sinh p. cosh ^ vr 



(7) 



27r 

The interaction among the elementary particles is attractive (repulsive) below (above) 
the free fermion point P = -\/47r. In the attractive regime, the elementary excitations 
form bound states (breathers). Our subsequent analysis can be easily extended to include 
breathers. For the sake of simplicity, however, we will not go into the details of the 
breather sector in this note. 

2. The basic physical ingredients entering the computation of form factors are conve- 
niently illustrated on the example of the two-fermion matrix element^ 

4(^^i -^2) = {O\m\A{0i)A,{e2)) , z,j = 1,2. (9) 

They satisfy the unitarity and crossing relations B 



Fm = E <(^)F*(-e) , 

fc,/=l,2 

F*(^ + 227r) = e2-^*.''F,f(-^) , (10) 



^The on-shell energy and momentum of the particles are parameterised as {p''\p^) 
(m cosh 9, m sinh 9) . 

^We denote by <i> a generic scalar (under Lorentz transformations) operator. 
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where the factor e^*'^'''*-* in the crossing equation takes into account a possible non-locahty 
among the operator $ and the particle 0, 0, |13| . For the operators we are concerned 



with in the AT model we just need to recall that Ai is associated to the fermion ijji. 
The latter is non-local with respect to the operators aj and fij with non-locality index 
1aj,i = 1fij,i = i^i is instead local with respect to itself and ej = ipjipj. For the 

operators defined as products, $ = Ilj'^i) the non-locality index reads 7$,j = 
These considerations provide the fourth column of the Table. 

It is obvious from the form of the lattice Hamiltonian (|l|) that the AT model is invariant 
under the change of sign of ai and/or a2 {Z2 x Z2-symmetry) , and under the exchange 
1^2 (exchange symmetry). The fermions ipi, and then the particles Ai, are odd under 
the Z2 X Z2-symmetry. The behaviour of the different operators under this symmetry then 
determines the asymptotic states on which each operator has non-zero matrix elements. 
These states are given in the last column of the Table (modulo pairs of the type AjAj). 
The symmetry properties of the operators under the exchange symmetry are important 
to simplify the functional system (|T0|). 

The scalar operators $ which are not odd under any of the global symmetries of the 
model may have a nonvanishing vacuum expectation value ($). If in addition they are 
non-local with respect to the fermions, their two-particle form factors exhibit a pole at 
9 = m with residue 

- iRese=i^Ff^{e) = (1 - e^^^^--){^) . (11) 

The only additional singularities allowed on the physical strip Q < 9 < m are the poles 
associated to the bound states. The functional equations (|To|) , together with the required 
singularity structure, fix the matrix elements (^ up to polynomials of cosh 9. The asymp- 
totic arguments of Ref. 0] can be used to show that, for the relevant (x<i, < 2) operators 
listed in the Table, these polynomials are of degree zero. The non-zero two-fermion matrix 
elements are then uniquely determined to be (j, / = 1, 2) 

= irr-, ■ . , [o{d) + i-ly^M2^rr - 9)] , (12) 

2t,uj[i7rj smh^(6' — 27rj l J 

Ff, = c^-J^?^;^F,i9), (13) 



sinh ^(9 — in) 



F^^=C2{-iyF,{9) , (14) 

Fli =F^^ =c^F^{9) , (15) 

Fk = -Fk = C4 . , ^o(^) , (16) 

cosh ^(t^ — ZTT) 

2^"'^^osh^(^-i7r) ' ^ 
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-T^^^inh^(^^-m) • ^^^^ 



Here the c„ are normalisation constants, (/i) = {^i) = {^2)1 and the functions 



dx sinh (1 - I) X sin 



,2 Ox 



• 2 (in—9)x 

sm ■'^ ^ 



6 i r°° dx sinh f f 1 - , 
Fo(^^) = -.sinh^exp / - . ^ ^ . , (20) 

2 I Jo X smh ^ cosh | smh x 



satisfy the equations 



sinh —(9-\- in) 

Lj{d) = uj{-9) , Luie + 2zn) = - . , ■ M d - 2^vr) , (21) 

Sinn 2^1'^ — ^'^ ) 

Fo{9) = S{9)Foi-9) , Fo{9 + 2zn) = Fo{-9) . (22) 

For large values of \9\ they behave as lu{9) ~ exp (^j — ij ^-^ and ^0(6*) ~ exp + 1^ . 

The form factors of the operators e*"*^ are known for the Sine-Gordon model 0, ^, 
Of course, the two-fermion matrix elements for the operators S, S^, V and V* determined 
here coincide with those of the corresponding bosonic operators in Sine-Gordon. 

The notation we use in the first column of the Table for the AT operators refers to the 
high temperature phase of the model. The high and low temperature phases, however, are 
mapped into one another by the duality transformations ai ^ fii, ipi ^ ipi. In particular, 
this implies that ai and fit represent the same operator in the two different phases. In the 
form factor approach their relative normalisation is fixed by the factorisation condition 

i 

hm Fti\9) = , (23) 

wich determines the one-particle matrix element Fp = (0|(jj(0)|Aj). 

Although the generalisation of the functional system ( p!ol) to the matrix elements with 
more than two fermions in the asymptotic state is straightforward, its solution poses a 
non-trivial mathematical problem 0. For the reasons to be explained in a moment, we 
do not need to go into this technical part here. 

3. In general, the knowledge of the ?7,-particle form factors (0|$(0)|n) allows expressing 
the correlation functions in the form of the spectral series 

00 

($i(x)$2(0)) = ^(0|<l>i(0)|n)(n|$2(0)|0)e-^"l^l , (24) 

n=0 

where En denotes the total energy of the n-particle asymptotic state. Although the exact 
resummation of the spectral series is beyond the present technical possibilities, it is known 
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that partial sums provide very good approximations (see e.g. Ref. for a more detailed 
discussion and additional references on this subject). Remarkably, this is already the case 
for the "two-particle approximation", which corresponds to neglecting the contribution 
of all the states with n > 2 in the spectral series (p^ . The fermionic matrix elements 
computed in this note are what is needed to implement the two-particle approximation 
in the breatherless region 47r < < Svr. As a quantitative check of the convergence of 
the spectral series, we evaluate in this approximation the central charge and the scaling 
dimension through the sum rules ll6|, H 
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(fx\x\^Q{x)Q{0))c, (25) 



X, = J d\ (e(x)^(O)), , (26) 

where (■ ■ ■)c denotes connected correlators and G(x) is the trace of the stress-energy tensor. 
The latter is proportional to the energy operator S{x) and its two-fermion form factors are 
given by ([13|) with the normalisation fixed byQ F^-{m) = 27rm^. Numerical integration 
shows that the two-fermion contribution to the central charge decreases monotonically 
from 1 at = 4tt to 0.987 at = Svr. For the two-fermion contribution to x^ one finds 
a monotonic decrease from 0.125 to 0.110 over the same range. The circumstance that 
this computation gives the exact answer at = Ait simply follows from the fact that 
<d ^ S = ipiipi + ip2'4'2 couples only to the two-particle states at the free fermion point. 
The fact that the approximation is more accurate for c than for away from the free 
point is due to the factor |xp in the integral (P^) which suppresses the contribution of the 
short distances. This is the region where our error localises when we truncate the large 



distance expansion (pif). 



Acknowledgements. I thank V.A. Fateev for discussions. 



^Our normalisation for the asymptotic states is defined by {A{9)\A{6')) — 2tt5{6 — 0'). 
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